Pyrazinoic acid, the active form of the anti-tuberculosis drug pyrazinamide, and aromatic carboxylic acid analogs are protonophores.

Pyrazinoic acid is the active form of pyrazinamide, a first-line antibiotic used to treat Mycobacterium tuberculosis infections. However, the mechanism of action of pyrazinoic acid remains a subject of debate, and alternatives to pyrazinamide in cases of resistance are not available. The work presented here demonstrates that pyrazinoic acid and known protonophores including salicylic acid, benzoic acid, and carbonyl cyanide m-chlorophenyl hydrazone all exhibit pH-dependent inhibition of mycobacterial growth activity over a physiologically relevant range of pH values. Other anti-tubercular drugs, including rifampin, isoniazid, bedaquiline, and p-aminosalicylic acid, do not exhibit similar pH-dependent growth-inhibitory activities. The growth inhibition curves of pyrazinoic, salicylic, benzoic, and picolinic acids, as well as carbonyl cyanide m-chlorophenyl hydrazone, all fit a quantitative structure-activity relationship (QSAR) derived from acid-base equilibria with R2 values > 0.95. The QSAR model indicates that growth inhibition relies solely on the concentration of the protonated forms of these weak acids (rather than the deprotonated forms). Moreover, pyrazinoic acid, salicylic acid, and carbonyl cyanide m-chlorophenyl hydrazone all caused acidification of the mycobacterial cytoplasm at concentrations that inhibit bacterial growth. Thus, it is concluded that pyrazinoic acid acts as an uncoupler of oxidative phosphorylation and that disruption of proton motive force is the primary mechanism of action of pyrazinoic acid rather than the inhibition of a classic enzyme activity.

Drug distribution and efficacy of the DprE1 inhibitor BTZ-043 in the C3HeB/FeJ mouse tuberculosis model.

BTZ-043, a suicide inhibitor of the Mycobacterium tuberculosis cell wall synthesis decaprenylphosphoryl-beta-D-ribose 2' epimerase, is under clinical development as a potential new anti-tuberculosis agent. BTZ-043 is potent and bactericidal in vitro but has limited activity against non-growing bacilli in rabbit caseum. To better understand its behavior in vivo, BTZ-043 was evaluated for efficacy and spatial drug distribution as a single agent in the C3HeB/FeJ mouse model presenting with caseous necrotic pulmonary lesions upon Mycobacterium tuberculosis infection. BTZ-043 promoted significant reductions in lung and spleen bacterial burdens in the C3HeB/FeJ mouse model after 2 months of therapy. BTZ-043 penetrates cellular and necrotic lesions and was retained at levels above the serum-shifted minimal inhibitory concentration in caseum. The calculated rate of kill was found to be highest and dose-dependent during the second month of treatment. BTZ-043 treatment was associated with improved histology scores of pulmonary lesions, especially compared to control mice, which experienced advanced fulminant neutrophilic alveolitis in the absence of treatment. These positive treatment responses to BTZ-043 monotherapy in a mouse model of advanced pulmonary disease can be attributed to favorable distribution in tissues and lesions, retention in the caseum, and its high potency and bactericidal nature at drug concentrations achieved in necrotic lesions.

Standardized RS Ratio Metrics To Assess Tuberculosis Antimicrobial Efficacy and Potency.

The sigmoid Emax model was used to describe the rRNA synthesis ratio (RS ratio) response of Mycobacterium tuberculosis to antimicrobial concentration. RS-Emax measures the maximal ability of a drug to inhibit the RS ratio and can be used to rank-order drugs based on their RS ratio effect. RS-EC90 is the concentration needed to achieve 90% of the RS-Emax, which may guide dose selection to achieve a maximal RS ratio effect in vivo.

Pharmacodynamics and Bactericidal Activity of Combination Regimens in Pulmonary Tuberculosis: Application to Bedaquiline-Pretomanid-Pyrazinamide.

A critical barrier to codevelopment of tuberculosis (TB) regimens is a limited ability to identify optimal drug and dose combinations in early-phase clinical testing. While pharmacokinetic-pharmacodynamic (PKPD) target attainment is the primary tool for exposure-response optimization of TB drugs, the PD target is a static index that does not distinguish individual drug contributions to the efficacy of a multidrug combination. A PKPD model of bedaquiline-pretomanid-pyrazinamide (BPaZ) for the treatment of pulmonary TB was developed as part of a dynamic exposure-response approach to regimen development. The model describes a time course relationship between the drug concentrations in plasma and their individual as well as their combined effect on sputum bacillary load assessed by solid culture CFU counts and liquid culture time to positivity (TTP). The model parameters were estimated using data from the phase 2A studies NC-001-(J-M-Pa-Z) and NC-003-(C-J-Pa-Z). The results included a characterization of BPaZ activity as the most and least sensitive to changes in pyrazinamide and bedaquiline exposures, respectively, with antagonistic activity of BPa compensated by synergistic activity of BZ and PaZ. Simulations of the NC-003 study population with once-daily bedaquiline at 200 mg, pretomanid at 200 mg, and pyrazinamide at 1,500 mg showed BPaZ would require 3 months to attain liquid culture negativity in 90% of participants. These results for BPaZ were intended to be an example application with the general approach aimed at entirely novel drug combinations from a growing pipeline of new and repurposed TB drugs.

Investigating the free-roaming dog population and gastrointestinal parasite diversity in Tulúm, México.

Free-roaming domestic dogs (Canis lupus familiaris) pose major conservation and public health risks worldwide. To better understand the threat of domestic dogs to wildlife and people and add to the growing literature on free-roaming dog ecology, a study was conducted to estimate the dog population in Tulúm, México. A modified mark-recapture technique and program MARK were used to obtain dog population estimates along six different transects dividing the city. Population estimates ranged from 19.75 dogs in one transect to 101.841 dogs in another, with 150 total dogs identified throughout the study and an estimated minimum population density of 48.57 dogs/km2. Fecal samples were also opportunistically collected for parasite identification through fecal flotation analysis using the McMaster technique. Out of 25 samples collected, 19 tested positive for gastrointestinal parasites with the most common species found being Ancylostoma caninum, followed by Toxocara canis, Dipylidium caninum, and Cystoisospora spp. Parasite loads ranged from 50 to 10,700 ova per gram of feces. The large population of free-roaming dogs and the prevalence of three zoonotic parasites highlight the importance of understanding free-roaming dog ecology and educating the public on the health risks free-roaming dogs pose. Los perros callejeros (Canis lupus familiaris) representan un gran riesgo para la conservación de animales y la salud pública mundialmente. Para comprender mejor la amenaza que significan los perros domésticos para la fauna silvestre y los humanos y aportar a la creciente bibliografía sobre la ecología de los perros callejeros, se realizó una investigación para estimar la población de los perros en Tulúm, México. Se utilizó una técnica modificada de marcado y recaptura junto con el programa MARK para estimar la población canina en seis transectos de la ciudad. Los estimados varían desde 19.75 perros en un transecto hasta 101,841 en otro, con un total de 150 perros identificados en el transcurso de la investigación y una densidad mínima estimada de 48,57 perros/km2. Además, se hizo una recolección oportunista de muestras de heces para la identificación de parásitos por medio del análisis de flotacíon fecal, con el método McMaster. De las 25 muestras recolectadas, 19 resultaron positivas para parásitos gastrointestinales, de las cuales las especies más comunes fueron Ancylostomoa caninum, seguida por Toxocara canis, Dipylidium caninum, y Cystoisospora spp. Las cargas parasitarias variaron desde 50 hasta 10.700 óvulos por gramo de heces. La alta población de perros callejeros y la prevalencia de tres enfermedades zoonóticas resaltan la importancia de entender la ecología de los perros callejeros y educar al público sobre los riesgos que significan los perros callejeros para la salud.

Pharmacodynamics and Bactericidal Activity of Bedaquiline in Pulmonary Tuberculosis.

Bedaquiline is a diarylquinoline antimycobacterial drug and a key component of several regimens in clinical development for the treatment of tuberculosis (TB) but with ongoing phase 3 trials that include assessment of simplified dosing. A pharmacokinetic-pharmacodynamic model of bedaquiline Mycobacterium tuberculosis-killing kinetics in adults with pulmonary TB was developed to inform dose selection of bedaquiline-containing regimens. The model parameters were estimated with data from the 14-day early bactericidal activity (EBA) study TMC207-CL001 conducted in Cape Town, South Africa. The study included 60 adult males and females with drug-susceptible pulmonary TB, who were administered bedaquiline with loading doses on the first 2 days followed by once-daily 100 mg, 200 mg, 300 mg, or 400 mg. The modeling results included expected values (means ± standard deviations [SDs]) for a maximum drug kill rate constant equal to 0.23 ± 0.03 log10 CFU/mL sputum/day, a half-maximum effective plasma concentration equal to 1.6 ± 0.3 mg/L, and an average time to onset of activity equal to 40 ± 7 h. Model simulations showed that once-daily 200 mg, 300 mg, and 400 mg (without loading doses) attained 40%, 50%, and 60%, respectively, of an expected maximum 14-day EBA equal to 0.18 log10 CFU/mL/day, or 10 h/day assessed by liquid culture time to positivity (TTP). Additional simulations illustrated efficacy outcomes during 8 weeks of treatment with the recommended and alternative dosages. The results demonstrate a general mathematical and statistical approach to the analysis of EBA studies with broad application to TB regimen development.

Pretomanid dose selection for pulmonary tuberculosis: An application of multi-objective optimization to dosage regimen design.

Clinical development of combination chemotherapies for tuberculosis (TB) is complicated by partial or restricted phase II dose-finding. Barriers include a propensity for drug resistance with monotherapy, practical limits on numbers of treatment arms for component dose combinations, and limited application of current dose selection methods to multidrug regimens. A multi-objective optimization approach to dose selection was developed as a conceptual and computational framework for currently evolving approaches to clinical testing of novel TB regimens. Pharmacokinetic-pharmacodynamic (PK-PD) modeling was combined with an evolutionary algorithm to identify dosage regimens that yield optimal trade-offs between multiple conflicting therapeutic objectives. The phase IIa studies for pretomanid, a newly approved nitroimidazole for specific cases of highly drug-resistant pulmonary TB, were used to demonstrate the approach with Pareto optimized dosing that best minimized sputum bacillary load and the probability of drug-related adverse events. Results include a population-typical characterization of the recommended 200 mg once daily dosage, the optimality of time-dependent dosing, examples of individualized therapy, and the determination of optimal loading doses. The approach generalizes conventional PK-PD target attainment to a design problem that scales to drug combinations, and provides a benefit-risk context for clinical testing of complex drug regimens.

Modeling and Simulation of Pretomanid Pharmacodynamics in Pulmonary Tuberculosis Patients.

Pretomanid (PA-824) is a nitroimidazole in clinical testing for the treatment of tuberculosis. A population pharmacodynamic model for pretomanid was developed using a Bayesian analysis of efficacy data from two early bactericidal activity (EBA) studies, PA-824-CL-007 and PA-824-CL-010, conducted in Cape Town, South Africa. The two studies included 122 adult male and female participants with newly diagnosed pulmonary tuberculosis who received once daily oral pretomanid doses of either 50, 100, 150, 200, 600, 1,000, or 1,200 mg for 14 days. The structural model described capacity-limited growth and saturable drug-induced bacterial killing with separate rate equations for sputum solid culture colony forming unit (CFU) counts and liquid culture time to positivity (TTP) that were linked through a time constant. The posterior population geometric means and interindividual variability percent coefficients of variation were, respectively; 0.152±0.013 log10 CFU/mL sputum/day and 54%±6% for the maximum kill rate constant, 20.4±1.0 h and 20.8%±0.1% for the time constant of proportionality between the CFU and TTP rate equations, and 770±140 ng/mL and 48%±17% for the pretomanid half-maximum effect plasma concentration. Model simulations showed once daily pretomanid at 100 mg, 200 mg, and 300 mg, attained 58%, 73%, and 80%, respectively, of an expected maximum 14-day EBA of 0.136 log10CFU/mL sputum/day. These results establish a pretomanid exposure-efficacy relationship with dual outcomes for CFU counts and TTP, and with potential applications to dose optimization of pretomanid-containing regimens.

Modeling and Simulation of Pretomanid Pharmacokinetics in Pulmonary Tuberculosis Patients.

Pretomanid is a nitroimidazole antibiotic in late-phase clinical testing as a component of several novel antituberculosis (anti-TB) regimens. A population pharmacokinetic model for pretomanid was constructed using a Bayesian analysis of data from two phase 2 studies, PA-824-CL-007 and PA-824-CL-010, conducted with adult (median age, 27 years) patients in Cape Town, South Africa, with newly diagnosed pulmonary TB. Combined, these studies included 63 males and 59 females administered once-daily oral pretomanid doses of 50, 100, 150, 200, 600, 1,000, or 1,200 mg for 14 days. The observed pretomanid plasma concentration-time profiles for all tested doses were described by a one-compartment model with first-order absorption and elimination and a sigmoidal bioavailability dependent on dose, time, and the predose fed state. Allometric scaling with body weight (normalized to 70 kg) was used for volume of distribution and clearance, with the scaling exponents equal to 1 and 3/4, respectively. The posterior population geometric means for the clearance and volume of distribution allometric constants were 4.8 ± 0.2 liters/h and 130 ± 5 liters, respectively, and the posterior population geometric mean for the half-maximum-effect dose for the reduction of bioavailability was 450 ± 50 mg. Interindividual variability, described by the percent coefficient of variation, was 32% ± 3% for clearance, 17% ± 4% for the volume of distribution, and 74% ± 9% for the half-maximum-effect dose. This model provides a dose-exposure relationship for pretomanid in adult TB patients with potential applications to dose selection in individuals and to further clinical testing of novel pretomanid-containing anti-TB regimens.

Bedaquiline and Pyrazinamide Treatment Responses Are Affected by Pulmonary Lesion Heterogeneity in Mycobacterium tuberculosis Infected C3HeB/FeJ Mice.

BALB/c and Swiss mice are routinely used to validate the effectiveness of tuberculosis drug regimens, although these mouse strains fail to develop human-like pulmonary granulomas exhibiting caseous necrosis. Microenvironmental conditions within human granulomas may negatively impact drug efficacy, and this may not be reflected in non-necrotizing lesions found within conventional mouse models. The C3HeB/FeJ mouse model has been increasingly utilized as it develops hypoxic, caseous necrotic granulomas which may more closely mimic the pathophysiological conditions found within human pulmonary granulomas. Here, we examined the treatment response of BALB/c and C3HeB/FeJ mice to bedaquiline (BDQ) and pyrazinamide (PZA) administered singly and in combination. BALB/c mice consistently displayed a highly uniform treatment response to both drugs, while C3HeB/FeJ mice displayed a bimodal response composed of responsive and less-responsive mice. Plasma pharmacokinetic analysis of dissected lesions from BALB/c and C3HeB/FeJ mice revealed that PZA penetrated lesion types from both mouse strains with similar efficiency. However, the pH of the necrotic caseum of C3HeB/FeJ granulomas was determined to be 7.5, which is in the range where PZA is essentially ineffective under standard laboratory in vitro growth conditions. BDQ preferentially accumulated within the highly cellular regions in the lungs of both mouse strains, although it was present at reduced but still biologically relevant concentrations within the central caseum when dosed at 25 mg/kg. The differential treatment response which resulted from the heterogeneous pulmonary pathology in the C3HeB/FeJ mouse model revealed several factors which may impact treatment efficacy, and could be further evaluated in clinical trials.

Computational pharmacokinetics/pharmacodynamics of rifampin in a mouse tuberculosis infection model.

One critical approach to preclinical evaluation of anti-tuberculosis (anti-TB) drugs is the study of correlations between drug exposure and efficacy in animal TB infection models. While such pharmacokinetic/pharmacodynamic (PK/PD) studies are useful for the identification of optimal clinical dosing regimens, they are resource intensive and are not routinely performed. A mathematical model capable of simulating the PK/PD properties of drug therapy for experimental TB offers a way to mitigate some of the practical obstacles to determining the PK/PD index that best correlates with efficacy. Here, we present a preliminary physiologically based PK/PD model of rifampin therapy in a mouse TB infection model. The computational framework integrates whole-body rifampin PKs, cell population dynamics for the host immune response to Mycobacterium tuberculosis infection, drug-bacteria interactions, and a Bayesian method for parameter estimation. As an initial application, we calibrated the model to a set of available rifampin PK/PD data and simulated a separate dose fractionation experiment for bacterial killing kinetics in the lungs of TB-infected mice. The simulation results qualitatively agreed with the experimentally observed PK/PD correlations, including the identification of area under the concentration-time curve as best correlating with efficacy. This single-drug framework is aimed toward extension to multiple anti-TB drugs in order to facilitate development of optimal combination regimens.

Computational pharmacology of rifampin in mice: an application to dose optimization with conflicting objectives in tuberculosis treatment.

Dose selection for rifampin in the treatment of active pulmonary tuberculosis (TB) illustrates some of the challenges for dose optimization within multidrug therapies. Rifampin-based anti-TB regimens are often combined with antiretroviral therapies to treat human immunodeficiency virus (HIV) coinfection. The potent cytochrome P450 (CYP) enzyme inducing properties of rifampin give rise to significant drug-drug interactions, the minimization of which by limiting the dose, conflicts with the maximization of bacterial killing by increasing the dose. Such multiple and conflicting objectives lead to a set of trade-off optimal solutions for dose optimization rather than a single best solution. Here, we combine pharmacokinetic/pharmacodynamic (PK/PD) modeling with multiobjective optimization to quantitatively explore trade-offs between therapeutic and adverse effects of optimal dosing for the example of rifampin in TB-infected mice. The PK/PD model describes rifampin concentrations in plasma and liver following oral administration together with hepatic CYP enzyme induction and bacterial killing kinetics. We include optimization objectives descriptive of antimicrobial efficacy, CYP-mediated drug-drug interactions, and drug exposure-dependent toxicity. Results show non-conventional dosing scenarios that allow for increased efficacy relative to uniform dosing without increasing drug-drug interactions. Additionally, we find currently employed dosages for rifampin to be nearly optimal with respect to trade-offs between efficacy and toxicity. While limited by the accuracy and applicability of the PK/PD model, these results provide an avenue for experimental investigation of complex dose optimization problems. This method can be extended to include additional drugs and optimization objectives, and may provide a useful tool for individualized medicine.

A physiologically based pharmacokinetic model of rifampin in mice.

One problem associated with regimen-based development of antituberculosis (anti-TB) drugs is the difficulty of a systematic and thorough in vivo evaluation of the large number of possible regimens that arise from consideration of multiple drugs tested together. A mathematical model capable of simulating the pharmacokinetics and pharmacodynamics of experimental combination chemotherapy of TB offers a way to mitigate this problem by extending the use of available data to investigate regimens that are not initially tested. In order to increase the available mathematical tools needed to support such a model for preclinical anti-TB drug development, we constructed a preliminary whole-body physiologically based pharmacokinetic (PBPK) model of rifampin in mice, using data from the literature. Interindividual variability was approximated using Monte Carlo (MC) simulation with assigned probability distributions for the model parameters. An MC sensitivity analysis was also performed to determine correlations between model parameters and plasma concentration to inform future model development. Model predictions for rifampin concentrations in plasma, liver, kidneys, and lungs, following oral administration, were generally in agreement with published experimental data from multiple studies. Sensitive model parameters included those descriptive of oral absorption, total clearance, and partitioning of rifampin between blood and muscle. This PBPK model can serve as a starting point for the integration of rifampin pharmacokinetics in mice into a larger mathematical framework, including the immune response to Mycobacterium tuberculosis infection, and pharmacokinetic models for other anti-TB drugs.

DoseSim: a tool for pharmacokinetic/pharmacodynamic analysis and dose reconstruction.

UNLABELLED: Assessing and improving the safety of chemicals and the efficacy of drugs depends on an understanding of the biodistribution, clearance and biological effects of the chemical(s) of interest. A promising methodology for the prediction of these phenomena is physiologically based pharmacokinetic/pharmacodynamic modeling, which centers on the prediction of chemical absorption, distribution, metabolism and excretion (pharmacokinetics) and the biological effects (pharmacodynamics) of the chemical on the organism. Strengths of this methodology include modeling across multiple scales of biological organization and facilitate the extrapolation of results across routes of exposure, dosing levels and species. It is also useful as the foundation for tools to (i) predict biomarker levels (concentrations of chemical species found in the body that indicate exposure to a foreign chemical), given a chemical dose or exposure; (ii) reconstruct a dose, given the levels of relevant biomarkers; and (iii) estimate population variability. Despite the importance and promise of physiologically based pharmacokinetic /pharmacodynamics-based approaches to forward and reverse dosimetry, there is currently a lack of user-friendly, freely available implementations that are accessible and useful to a broad range of users. DoseSim was developed to begin to fill this gap. AVAILABILITY: The application is available under the GNU General Public License from http://scb.colostate.edu/dosesim.html.

Computational toxicology of chloroform: reverse dosimetry using Bayesian inference, Markov chain Monte Carlo simulation, and human biomonitoring data.

BACKGROUND: One problem of interpreting population-based biomonitoring data is the reconstruction of corresponding external exposure in cases where no such data are available. OBJECTIVES: We demonstrate the use of a computational framework that integrates physiologically based pharmacokinetic (PBPK) modeling, Bayesian inference, and Markov chain Monte Carlo simulation to obtain a population estimate of environmental chloroform source concentrations consistent with human biomonitoring data. The biomonitoring data consist of chloroform blood concentrations measured as part of the Third National Health and Nutrition Examination Survey (NHANES III), and for which no corresponding exposure data were collected. METHODS: We used a combined PBPK and shower exposure model to consider several routes and sources of exposure: ingestion of tap water, inhalation of ambient household air, and inhalation and dermal absorption while showering. We determined posterior distributions for chloroform concentration in tap water and ambient household air using U.S. Environmental Protection Agency Total Exposure Assessment Methodology (TEAM) data as prior distributions for the Bayesian analysis. RESULTS: Posterior distributions for exposure indicate that 95% of the population represented by the NHANES III data had likely chloroform exposures < or = 67 microg/L [corrected] in tap water and < or = 0.02 microg/L in ambient household air. CONCLUSIONS: Our results demonstrate the application of computer simulation to aid in the interpretation of human biomonitoring data in the context of the exposure-health evaluation-risk assessment continuum. These results should be considered as a demonstration of the method and can be improved with the addition of more detailed data.